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Abstract 

An evaluation of several solar concentrator-based systems for producing oxygen from lunar regolith was 
performed. The systems utilize a solar concentrator mirror to provide thermal energy for the oxygen production 
process. Thermal energy to power a Stirling heat engine and photovoltaics are compared for the production of 
electricity. The electricity produced is utilized to operate the equipment needed in the oxygen production process. The 
initial oxygen production method utilized in the analysis is hydrogen reduction of ilmenite. Utilizing this method of 
oxygen production a baseline system design was produced. This baseline system had an oxygen production rate of 
0.6 kg/hr with a concentrator mirror size of 5 m. Variations were performed on the baseline design to show how 
changes in the system size and process (rate) affected the oxygen production rate. An evaluation of the power 
requirements for a carbothermal lunar regolith reduction reactor has also been conducted. The reactor had a total power 
requirement between 8,320 to 9,961 W when producing 1000 kg/year of oxygen. The solar concentrator used to 
provide the thermal power (over 82 percent of the total energy requirement) would have a diameter of less than 4 m. 


Introduction 

The ability to produce oxygen is a key aspect of the successful human exploration and colonization of the lunar 
surface. Oxygen is required as a propellant reactant and is also needed to sustain a breathable atmosphere within 
habitation modules. Transporting oxygen from the Earth’s surface is very costly and requires a sufficient stockpile to be 
housed on the lunar surface as a precaution in case there is a disruption with the supply from Earth. By producing the 
oxygen from material available on the lunar surface, the dependence on Earth resupply and its cost can be eliminated. 

In addition to the cost savings, utilizing available resources to produce oxygen is the first step toward a viable self- 
sustaining human presence on the lunar surface. 

Oxygen is the most prevalent element on the lunar surface, about 45 percent of lunar regolith (Ref. 1). However, it 
is chemically locked up in various compounds that make up the regolith. There are a number of reactions and processes 
that have been devised which can break apart the regolith compounds and liberate the oxygen (Ref. 2). One common 
aspect of all of these processes is the need for heat to be added to the regolith. The most efficient method of achieving 
this heat addition is by direct heating of the regolith with concentrated solar energy. 

This study examines two cases: a system that utilizes a solar concentrating mirror as a means of providing the 
thermal energy to heat the regolith and a Stirling heat engine to provide the electrical power to operate the processes 
involved in the production of oxygen. There is an inherent advantage in overall system efficiency and compact system 
design by providing direct thermal heating as well as power generation from a single concentrator mirror. The ultimate 
goal of this kind of compact arrangement would be the design of an autonomous, self-contained oxygen production 
system that can be deployed to the lunar surface and begin oxygen production without the need for human assembly. 
The other system employs a hybrid approach with photovoltaics to produce electrical energy. Greater design flexibility 
through choice of solar cell technology is gained through this more complex approach. 
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Nomenclature 


A s surface area 

c p ilmenite specific heat 

d t ilmenite maximum particle size 

du spacing between insulation layers 

d m concentrator mirror diameter 

d r reactor diameter 

d rs regolith transport belt distance 

g surface gravitational constant 

H r heat of reaction 

h r reactor height 

h rs reactor height above the surface 

I m mean solar intensity 

M. mass rate of ilmenite to the reactor 

M o rate of oxygen production 

m Q desired oxygen production per year 

M r mass rate of regolith acquisition 

• 

M w rate of water production 

rii number of layers of insulation 

n m number of melts being processed 

n r regolith particle size distribution 

P e electrolyzer power 

Pm hydrogen gas pressure 

Pu power to lift the ilmenite 

Pms regilith/ilmenite separator power 

P rs regolith scoop power 

P r i power to transport the regolith to the screen 

P se total required power 

P vs vibratory screener power 

Q hi heat required to raise the ilmenite to the reactor temperature 

Qi s heat loss to the surroundings 

Q r heat consumed in the reaction 

Q se heat required by the Stirling engine 

gtotai total heat required 

Pm hydrogen gas constant 

R m ilmenite mass to regolith mass ratio 

r 0 rb orbital radius of earth 

''orbm mean orbital radius of earth 

th ilmenite heating time 

Ti ilmenite temperature within the reactor 

T io ilmenite initial surface temperature 

T m regolith melt temperature 

t m time needed to melt regolith 

T r regolith initial temperature 

T s temperature of the surroundings 

t s total available sunlight time 

t t ilmenite time at reactor maximum temperature 

U m hydrogen gas velocity within the reactor 

• 

Vir reactor volume replaced each hour 

v m volume of a single melt 

y Q total oxygen yield 
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Z H 2 hydrogen gas compressibility factor 

a solar day angle 

e eccentricity of earth’s orbit 

8 i insulation layer emissivity 

z p porosity of the fluidized bed 

z r reactor surface emissivity 

\\ e electrolyzer efficiency 

r\ em electric motor efficiency 

r\ m concentrator mirror reflectance efficiency 

v\ r heat receiver efficiency 

r\ s efficiency of the separation process 

r\ se efficiency of the Stirling engine 

Hh 2 hydrogen gas viscosity 

Hr regolith coefficient of friction 

p j ilmenite density 

p ip ilmenite particle density 

Ph 2 hydrogen gas density 

pr regolith density 

a Boltzman constant 

cp ilmenite particle shape factor 


The Lunar Environment 

The Moon (shown in Fig. 1) (Ref. 3) is a dry, crater-filled land with no appreciable atmosphere. The 
temperature difference from the sunlit areas to shadowed areas can exceed 300 K. The surface is covered with a fine 
powdered regolith made from eons of bombardment by meteors, asteroids, and comets. The successful design and 
operation of a lunar oxygen production system within this harsh environment will require the system to leverage the 
advantages of the environment (such as abundant solar energy) and minimize the disadvantages (such as temperature 
extremes). From the solar availability to the amount of removable oxygen in the soil, the capabilities of this system 
will be based on how well it is adapted to these environmental conditions. 



Figure 1. — Image of the Moon. 
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There are two main types of terrain on the lunar surface, the highlands and the mare. The highlands are the 
older, brighter regions (as seen from Earth) of the lunar surface and are significantly cratered. The mare consist of 
darker regions which are basaltic lava flows. These areas can be seen in the topographic map of the lunar surface 
shown in Figure 2. The blue and purple regions are the lower-lying mare regions while the red, orange and yellow 
areas are the upper highland regions. The lunar regolith, which covers the surface of the Moon, consists of fragments 
of rocks, minerals and glass spherules formed when meteors or other bodies impacted the surface. The regolith 
thickness varies greatly with location on the lunar surface. In the mare the thickness can vary between 3 to 16 m 
whereas in the highlands it is at least 10 m thick (Ref. 1). The chemical or mineral composition of the regolith is 
similar to the underlying bedrock from which it was derived. There is little mixing of material between the highland 
and mare regions. The composition of the major materials (greater than 1 percent of the composition) in each of 
these regions is given in Table 1. 

The regolith is composed of fine grains of the mineral compounds given in Table 1. The grain size of the 
regolith varies but consists mostly of particles less than 1 mm in diameter. This is similar in consistency to silt or 
fine sand (Ref. 1). Average regolith density (p r ) is approximately 1000 kg/m 3 . Key compounds in the production of 
oxygen from lunar material by hydrogen reduction are the iron oxides and titanium dioxides. 

The highest concentrations of iron oxide and titanium dioxide (and ilmenite, FeTi0 3 ) occur in the low-lying 
mare regions. Within these mare regions of high iron and titanium content, ilmenite will be present in the loose 
regolith. About two thirds of the ilmenite is in glassy agglutinates and basaltic rock fragments much larger than the 
common regolith grain size. The average amount of small grain size ilmenite available in the mare regions is 
approximately 7 percent by weight (the ratio of mass of ilmenite to mass of regolith, R m , is 0.07) and 5 percent by 
volume (the ratio of volume of ilmenite to volume of regolith, R v , is 0.05) of the regolith (Ref. 4). Based on these 
percentages, the density of ilmenite (p t ) is approximately 1.4 times that of the regolith or 1400 kg/m 3 . The thermal 
conductivity of the fine-grained ilmenite is estimated to be 0.425 W/mK which is similar to that of sand (Ref. 5). 



Figure 2. — Lunar topographic map from Reference 3. 


TABLE 1.— LUNAR SURFACE MINERAL COMPOSITIONS 
(REFS. 1 AND 2) 


Compound 

Highland, 

Mare, 


wt% 

wt% 

Si0 2 

44.5 

41.0 

ai 2 0 3 

26.0 

12.8 

FeO 

5.77 

16.2 

CaO 

14.9 

12.4 

MgO 

8.05 

9.2 

Ti0 2 


7.3 
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Figure 3. — Variation in solar intensity at Earth orbital location. 


The specific heat of ilmenite is approximated by that of the mare regolith. The specific heat ( c p ) as a function of 
temperature ( T) is given by Equation (1) (Ref. 5). 

c p = -1848.5+1 047.41 log(7J) [J/kgK] (1) 

With the absence of an appreciable atmosphere, the solar intensity at the lunar surface is fairly constant. There 
is a slight variation throughout the year due to the Earth’s orbital eccentricity about the Sun. The mean solar 
intensity (7 m ) at Earth’s orbital location is 1353 W/m 2 (Ref. 6). The variation in Earth’s orbital radius (r orb ) from the 
mean orbital radius (r orbm , which has a value of 1 .496E8 km) is represented by the eccentricity (e) of Earth’s orbit 
that has a value of 0.017. The actual solar flux (or intensity, I) in W/m 2 for a specific day of the year is determined 
by Equations (2) and (3). This variation in solar intensity throughout the year is shown in Figure 3. 

7=4 (^orbm^orb 2 ) (2) 

rorb = (l-e 2 )/0+eCos(a)) (3) 

The day angle (a) used in Equation (3) is defined as 0° on January 4th (perihelion of Earth’s orbit) and 
increases by 0.98° per day. The solar output varies from approximately 1400 W/m 2 (on January 4th) to 1308 W/m 2 
(on July 9th). The variation of the solar intensity at the lunar surface due to the orbit of the Moon around the Earth is 
very slight. The distance between the Moon and Earth varies between 360,300 km (at perigee) to 405,500 (at 
apogee). The greatest effect on solar intensity produced by the Moon’s orbit about the Earth is a decrease of about 
0.56 percent (~6 W/m 2 ) in the solar intensity at the Lunar surface. This small change in intensity will vary as the 
Moon orbits the Earth or in other words throughout the lunar day. Therefore it was assumed to be negligible and the 
solar intensity at Earth orbit was utilized in the subsequent analysis. 



Oxygen Production by Hydrogen Reduction of Ilmenite — Process Overview 

There are a number of potential methods for extracting oxygen from the lunar soil (Ref. 2). The method 
examined in this section of our analytical study is the hydrogen reduction of ilmenite. This is one of the methods of 
oxygen production from Lunar regolith that is under active development. It is a fairly simple method that can be 
accomplished with present day materials. Because the reaction takes place at a relatively low temperature, below the 
melting point of the soil, the material compatibility and corrosiveness are not significant issues. The process utilizes 
hydrogen gas that will need to be supplied from Earth. However, the gas is recoverable in the oxygen generation 
process. Therefore after the initial supply only make-up gas will be needed to account for that which is not 
recovered or lost through leaking. Another advantage of this type of production is that it can utilize the surface 
regolith in areas that have a high iron and titanium content. This eliminates the need for the mining or processing of 
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rocks, these are both very power intensive. For an initial oxygen production capability, and as a means of 
demonstrating the oxygen production technology, being able to directly utilize regolith is a significant advantage. 

Hydrogen reduction of ilmenite is a two-step process represented by the reactions given in 
Equations (4) and (5) (Refs. 2 and 4). The first step in the process is the reduction of the ilmenite to iron, titanium 
dioxide, and water. The second step is the electrolysis of the water to form hydrogen and oxygen. 

2 FeTi0 3 + 2 H 2 ^ 2 Fe + 2 Ti0 2 + 2 H 2 0 (4) 

2 H 2 0 ■> 2 H 2 + 0 2 (5) 

Both the hydrogen and oxygen that are produced in the process can be captured. The oxygen will be stored for use 
and the hydrogen will be reused in the process. The reduction process is performed in a fluidized bed reactor and the 
electrolysis is performed with an electrolyzer. A diagram of the process is given in Figure 4. 

The process requires the ilmenite to remain at an elevated temperature (7)) of 1273K for approximately 1 hr. An 
additional benefit to the heating of the ilmenite is the liberation of hydrogen. Hydrogen is present in mature lunar 

regolith due to exposure to the solar wind over long periods of time. The finer the particle grains the larger the 

amount of hydrogen it will have trapped. Hydrogen content will range from 0.1 to 0.2 kg/m 3 of regolith. Since 
hydrogen content is dependent on grain size not material content it is reasonable to assume the concentration in the 
ilmenite is the same as that in the overall regolith. Therefore, as the ilmenite is heated the hydrogen bonded to it will 
be released. It is estimated that at 900 °C all of the hydrogen trapped with the ilmenite will be released. Released 
hydrogen will be utilized in the reduction of the ilmenite and will add to the overall hydrogen already in the system. 
It can be considered a source of makeup hydrogen to help offset any losses of hydrogen from the system. The 
sequence for the process illustrated in Figure 4 consists of the following: 


• Scooping the regolith and transporting it to the screen. 

• The larger regolith particles (greater than 1 mm in diameter) are separated using a vibrating screen. 

• The regolith is then transported to a magnetic separator where the ilmenite is separated from the regolith. 

• The ilmenite is transported to the fluidized bed reactor. 

• In the reactor the ilmenite is heated while hydrogen gas is circulated (it is maintained at 1273 K for 1 hr). 

• The water vapor output is passed through an electrolyzer where it is separated into hydrogen and oxygen gas. 

• The hydrogen gas is returned to the reactor and the oxygen is stored. 


concentrator 

mirror 


Regolith 

belt 

scoop — 



— Oxygen tank 


Hydrogen tank 
"" — Stirling heat engines 

Large ilmenite Depleted ilmenite 

particles separated >- Regolith after ( Fe and Ti0 2) 

ilmenite is removed 


Figure 4. — Diagram of hydrogen reduction of ilmenite oxygen production 
processing system. 
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Methodology for Sizing of Integrated Power System for Hydrogen Reduction 

An analysis was set up to evaluate the sizing of the system shown in Figure 4 for various rates of oxygen 
production (Ref. 7). 

The sizing was based on scaling the solar dynamic system to provide sufficient heat to run the ilmenite 
reduction reaction as well as provide electrical power for operating all of the other process involved with the oxygen 
production. The electrical power is produced by a series of Stirling heat engines. To perform the sizing analysis an 
energy balance between the heat provided by the concentrator mirror and the requirements of the various 
components both thermally and electrically had to be made. Each of the components’ thermal or electrical energy 
consumption has to be determined and scaled with the amount of regolith and or ilmenite being processed. Since the 
process involves moving the material around, the height and rate of movement is needed in order to determine the 
energy consumption. These values can be varied to determine what effect they have on the system sizing and 
production rate. 

The rate of ilmenite going into the reactor sets the rates for all of the other processes. This rate can be 
represented by the fraction of the reactor volume that is replenished each hour, v ir • This quantity will have a value 
between 0 and 1 . Establishing the replacement rate of the ilmenite in the reactor also sets the volume and mass flow 
requirements for all the other process in the system. From this, the mass rate at which the regolith must be acquired 
from the surface ( M r in kilograms per hour) can be determined. This rate, which is based on the total volume of 
ilmenite in the fluidized bed reactor, is given by Equation (6). The reactor is assumed to be cylindrical. The diameter 
(dy) and height (h r ) can be adjusted to increase or decrease the overall volume of ilmenite processed. 


M = 


p i V if .nd r h r 
4 R m 


( 6 ) 


The rate at which the regolith is taken from the surface sets the power required to scoop the regolith (P rs ) and 
transport it to the screen (P r ). These quantities are given in Equations (7) and (8), respectively. The power to scoop 
the regolith is dependent on the distance the belt scoops travel in the regolith (< d rs , assumed to be 1 m) and the 
coefficient of friction of the regolith (p r ). The regolith coefficient of friction was assumed to be 0.4, which is similar 
to that of sand. The power to lift the regolith is dependent on the height it is raised. This height is assumed to be the 
height of the reactor plus the height the reactor is above the surface (h rs , assumed to be 1 m). The electric motor 
efficiency ( x \ em ) to drive the scoop and lift the regolith is assumed to be 90 percent. 


Mrejfr 8 

3600n em 


( 7 ) 


(h r +h rs )M r g 
3600n em 


( 8 ) 


Once the regolith is picked up off of the surface, the first step is to screen out the larger particles; those greater 
then a specific size are removed. This screening is necessary for the efficient operation of the fluidized bed reactor. 
The maximum particle size (d s ) and subsequent screen opening that allows for efficient fluidization is calculated in 
Equation (9) (Ref. 4). This equation has to be solved iteratively for the maximum particle size. It factors in a number 
of parameters for the operation of the fluidized bed reactor. These include the porosity of the fluidized bed (e^, 
estimated to be 0.5) and the hydrogen gas velocity within the reactor ( U m > estimated to be 0.3 m/s). Details on the 
hydrogen gas properties and the ilmenite are also required. These include the hydrogen gas density (pn 2 , given in 
Equation (10)), hydrogen viscosity (p,H 2 , estimated to be 2.37E-5 kg/m-s), (Ref. 4) ilmenite particle shape factor (c|), 
estimated to be that of sand, 0.83) and the ilmenite particle density (p ip9 estimated to be 4790 kg/m 3 ). The ilmenite 
particle density is greater than that of the ilmenite in the soil. This is because the density of the ilmenite in the soil 
takes into account any spacing between the particles where as the particle density is that of the material that makes 
up each particle. 
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( 9 ) 


150p Hi(}-Zp)UH2 SP mip ip ~Phi) ^2 1.75p# 2 ^2 ^ 

<P 2 Z 3 pVH2 VH2 1 < p 4^2 

The hydrogen gas density is based on the ideal gas law (pressure (P H 2 ) in Pascal, temperature (7}) in degrees 
Kelvin and the gas constant for hydrogen (Rm), 4157 Nm/kgK) with a compressibility factor (Z H2 ) for the hydrogen 
gas (Ref. 8). 


Ph2 - 


Ph2 

Z U2 R m T i 


( 10 ) 


Z m = 0.99704 + 6.4149E-9P m 


( 11 ) 


The screening method selected for this analysis is a vibratory screen. A vibratory screening process is fairly 
simple. A screen with a specified mesh or opening size is vibrated as the regolith is poured onto it. The material that 
is less than the opening size will fall through to the next stage in the process whereas the material larger than the 
openings will be moved across the screen surface and discharged off one end of the screen and returned to the lunar 
surface. The power required to operate the vibratory screener (P vs ) is given by Equation (12) (Ref. 4). 


0.45 M r 

507.72*/, + 49.91^ 


( 12 ) 


The screening process will eliminate a portion of the regolith that was collected. The fraction of the regolith that 
will make it through the screen will depend on the desired maximum particle size determined by Equation (9). The 
particle size distribution (n r ) of the regolith in meters is given by Equation (13). This equation is a curve fit 
representation of data on the percent concentration of particle size within the regolith (Ref. 4). The particle size 
distribution along with the corresponding curve fit is plotted in Figure 5. 

d t -9.4008ZT-6 

n r =— (13) 

rf,+8.0419£-5 



Particle size 


Figure 5. — Particle size distribution of Lunar regolith. 
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The next step in the process is to separate out the ilmenite from the regolith. For a smaller scale oxygen 
production plant the use of a permanent or electro-magnet separator is preferred. This type of separator minimizes 
the power requirements compared to an induced magnetic roll separator or an electrostatic separator. This type of 
separator consists of a drum with alternating magnetic poles along its surface. Material is dropped over the drum as 
it rotates. Magnetic material will adhere to the surface where it is scrapped off and captured. Nonmagnetic material 
will drop out of the bottom of the separator (Ref. 9). Power consumed by the separator (P^) in Watts is given by 
Equation (14) (Ref. 4). 


Pms =0.19 6n r M r 


(14) 


The separated ilmenite has to now be lifted into the fluidized bed reactor. This distance is assumed to be the 
height of the reactor plus the distance the reactor is above the surface. The power required to lift the ilmenite (Pu) is 
given by Equation (15). Where the mass flow of ilmenite in kilograms per hour is given by Equation (16). It is 
assumed that the separation process is 90 percent efficient ( r \ s ). Therefore 90 percent of the mass flow given by 
Equation (15) will be ilmenite. The remainder will be other materials that will not contribute to the reaction although 
it will go through the complete heating process. 


p _ (K+KsWjg 

U 3600n em 


(15) 


M t = n r R m M r 


(16) 


Once the ilmenite enters the fluidized bed reactor it is heated to 1273 K and hydrogen gas is circulated through 
it. The energy to heat the soil comes from the solar concentrator. The concentrated solar radiation is focused onto a 
heat receiver. This heat receiver then conducts the heat into the reactor chamber where it is distributed throughout 
the chamber to heat the ilmenite. It is assumed that the ilmenite will heat uniformly once it is placed in the reactor. 
As the ilmenite passes from the top of the reactor to the bottom it is slowly heated. The selected rate, at which the 
volume of ilmenite in the reactor is replaced, given previously as v i , sets the amount of heat needed to heat the soil. 
The ilmenite must be brought up to temperature and remain at temperature for a specified period of time (t t9 in 
hours) in order for the reduction reaction (given by Eq. (4)) to take place. The time available to heat the ilmenite 
(f*, in hours) is given by Equation (17). From this equation it should be observed that there is an inherent limit on the 
rate at which the ilmenite is replaced and the duration at which it must be at temperature. This is because the time 
available to heat the ilmenite must be greater than zero. 

**=(4— t t ) where V ir t, < 1 (17) 

*l> 


The thermal power required to heat the ilmenite ( Q hi ) , given by Equation (18), is dependent on this available 
heating time, the volume of ilmenite being heated and the initial ilmenite properties. It is assumed that the ilmenite is 
initially at the illuminated surface temperature ( T io ) of 384 K. 


Qhi - 


P i c P (Ti ~ T jo )nd?h r (1 -t t Vjr) 

14,400f A 


(18) 


The ilmenite will lose heat to the surroundings (Qi s ) as it is being heated. To minimize this heat loss, the outside 
of the reactor vessel can be insulated with multilayer insulation. Because of the vacuum conditions at the lunar 
surface, radiation heat transfer is the only mechanism in which the reactor vessel will transfer heat to the 
surroundings. Therefore multiple layers of a low emissivity (e t ) material (such as polished nickel, = 0.17) that are 
closely layered will provide very good insulation. The greater the number of layers (nj) and the smaller the spacing 
between the layers (dn) 9 the less the heat loss to the surroundings will be. The heat loss to the surroundings (at a sink 
temperature T s = 270 K) from the insulated reactor vessel can be calculated from Equation (19) where the Stefan- 
Boltzman constant (a) has a value of 5.67E-8 W/m 2 K. If the reactor is constructed of a material different than that of 
the insulation it will have a different emissivity (e r . such as chromium e r = 0.14). 
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(19) 


Qis = 


^d r+ni d a )hMT? -T*) 


B „ 8, 


In addition to the heat required to bring the ilmenite to the desired temperature and the loss of heat to the 
surroundings, the hydrogen reduction reaction itself will absorb heat from the system. The amount of heat consumed 
( Q r ) due to the heat of reaction (/£.) is given by Equation (20). The heat of reaction for the hydrogen reduction of 
ilmenite is given as 294,000 J/kg (Ref. 4). 


Qr = 


H r T\ s Mi 

3600 


( 20 ) 


The last power-consuming component of the system is the electrolyzer. The electrolyzer power requirement will 
depend on the amount of water that needs to be electrolyzed. The rate of water production (m w ) is based on the ratio 

of the molecular weights of the reactant and product of interest, water. The rate of water production in kilograms per 
hour is given by Equation (21). Subsequently the rate of oxygen production (m q ) can then be determined, which is 
given by Equation (22). 


( 21 ) 

M a =M w ~ ( 22 ) 

The power required by the electrolyzer to separate the water into hydrogen and oxygen can be calculated based 
on the theoretical power needed to break apart water (3520 W/kg/hr), the efficiency of the electrolyzer (r\ e ) and the 
rate of water flow to the electrolyzer, see Equation (23). 

P e = 3520^- (23) 

n e 

The electrical power to operate the scoop and belt motors, run the magnetic separator and vibratory separator 
and operate the electrolyzer is provided by a number of Stirling heat engines. The required power output (P se ) of 
these engines is given by Equation (24) and the subsequent heat input to the system to run the engines is given by 
Equation (25). This heat input is based on the conversion efficiency of the Stirling engine (p je ) that is estimated to 
be 35 percent. 


P se - P rs + P rl + P vs + P ms + P il + P e (^ 4 ) 

Qse=~ ( 25 ) 

T\se 

The total thermal power output needed and the corresponding diameter of the solar concentrator mirror (d m ) can 
now be calculated. The total thermal power (Qtotai) output is given by Equation (26). The required solar concentrator 
diameter is given by Equation (27) and assumes it is circular in shape. The diameter is dependent on the mirror 
reflectance efficiency (r| m ) and the efficiency of the heat receiver (q r ). 

QxoXsi = Qse + Qr + Qls + Qhi ( 26 ) 


r 1 


( 27 ) 
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Results of Sizing of Integrated Power System for Hydrogen Reduction 

Utilizing the analysis outlined in the previous section, data was produced on the achievable oxygen production rate 
for various size systems. Initially a baseline case was produced that represents the initial design and operation point for 
the system. The baseline assumptions that were used to generate this case are given in Table 2. These assumptions 
represent the best estimates of the performance for the various components as well as an initial sizing that is within the 
scope of a demonstration oxygen production system. The system sizing and performance results for the baseline system 
are given in Table 3; a breakdown of the thermal and electrical power requirements is shown in Figure 6. 


TABLE 2.— BASELINE ANALYSIS ASSUMPTIONS 


Variable 

Value 

Fraction of the fluidized bed reactor that is replenished each hour ( V ir ) 

0.1 (1/hr) 

Reactor diameter (d r ) 

0.5 m 

Reactor height (h r ) 

0.5 m 

Operational date, day angle (a) 

July 9th, 180° 

Distance the scoop is in the Regolith (d rs ) 

1 m 

Height the fluidized bed reactor is above the surface ( h 

1 m 

Porosity of the fluidized bed (e^,) 

0.5 

Hydrogen velocity through the fluidized bed ( U m ) 

0.3 m/s 

Hydrogen pressure within the fluidized bed (Pm) 

1E6 Pa (10 atm) 

Ilmenite hydrogen reduction reaction temperature (T t ) 

1273K 

Time ilmenite is held at reaction temperature (t t ) 

1 hr 

Number of layers of MLI insulation («/) 

50 

Sink temperature of surroundings (T s ) 

270K 

Electric motor efficiency (r| em ) 

0.90 

Ilmenite separation efficiency (r\ s ) 

0.90 

Electrolyzer efficiency 

0.72 

Stirling engine efficiency 

0.35 

Concentrator mirror efficiency 

0.90 

Heat receiver efficiency 

0.90 


TABLE 3.— BASELINE SIZING AND PERFORMANCE RESULTS 


Quantity 

Value 

Available solar intensity 

1301 W/m 2 

Total thermal power required 

15,241 W 

Total electrical power required 

3,272 W 

Concentrator mirror area/diameter 

14.38 m 2 /4.28 m 

Reactor volume 

0.098 m 3 

Maximum particle size 

0.842 mm 

Total mass of hydrogen required 

0.73 kg 

Water production rate 

0.653 kg/hr 

Oxygen production rate 

0.58 kg/hr 
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Heat loss to surroundings 
■ llmenite heating 
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Operate electrolyzer 
■ Process Regolith 
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Figure 6. — Breakdown of thermal power (a) and electrical power (b) requirements. 



Ilminate volume, m 3 llminate volume, m 3 


Figure 7. — Oxygen production rate and mirror diameter Figure 8. — Thermal and electrical output requirements 

for various volume reactors. for various volume reactors. 

From this initial case, variations were made to some of the variables to determine what impact they have on the 
system sizing and output. The container volume was varied from near zero to 2.5 m 3 . The remaining variables, given 
in Table 2, were not changed. The corresponding oxygen production rate and mirror diameter are plotted in Figure 7 
and the required thermal and electrical output is plotted in Figure 8. From these figures it can be seen that there is a 
linear increase in oxygen production with the increase in reactor volume. This should be expected since the amount 
of oxygen is directly proportional to the volume of ilmenite being reduced. The thermal and electrical power 
requirements also increase linearly with the reactor volume. 

Another means of changing the rate of oxygen production is to vary the rate at which the ilmenite is replaced in 
the reactor. For the baseline case this rate was set at 10 percent of the container volume per hour. By varying this 
rate the amount of oxygen produced will change as well as the required thermal power and correspondingly mirror 
size. The oxygen production rate and the required mirror diameter as a function of the fraction of the ilmenite in the 
reactor replaced each hour is shown in Figure 9. Figure 10 shows the thermal and electrical power required as a 
function of the fraction of ilmenite replaced in the reactor each hour. 

Using the data produced a direct comparison between the oxygen production rate and the mirror diameter can 
be made. This comparison is plotted in Figure 1 1 for both the variable volume and the variable rate cases. From this 
figure it can be seen that for a given oxygen production rate the required mirror diameter is basically the same 
between the variable volume or rate cases. Therefore higher production rates can be achieved either through 
increasing the volume of the reactor or by increasing the rate at which the ilmenite is processed through the reactor. 
In either case, for a given rate of oxygen production the effect on the thermal and power system sizing is the same. 
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Faction of ilminate replaced per hour 


Figure 9. — Oxygen production rate and mirror diameter 
for different rates of ilmenite replacement. 



Faction of ilminate replaced per hour 

Figure 10. — Thermal and electrical output requirements 
for different rates of ilmenite replacement. 



0 1 2 3 4 5 6 

Oxygen production rate, kg/hr 

Figure 1 1 . — Mirror size as a function of oxygen 
production rate for both variable volume and rate. 


Oxygen Production by Carbothermal Processing 

The Carbothermal Reduction process utilizes carbon to bond with the oxygen contained in the oxides within the 
regolith. Each of the compounds in Table 1 (materials less than 1 percent were excluded) reacts under high 
temperatures (1900 °C) in the presence of a carbon rich environment (in this case, methane, CH4) (Ref. 10). 

Methane is introduced to the metal oxides at the elevated temperature, and extracts oxygen (using the carbon) to 
leave behind the metal (Eq. (28)). The CO and H 2 are then transferred to a methanation reactor where additional H 2 
is added. These gases react to form CH 4 (which would be reused in the first reaction) and H 2 0 (Eq. (29)). The last 
step would involve electrolyzing the H 2 0 to form H 2 (which would be reused in the methanation reaction) and 0 2 
(which would be stored) (Eq. (30)). 


MO x + x CH 4 -> M + x CO + 2x H 2 

(28) 

x CO x + 3x H 2 -> x CH 4 + x H 2 0 

(29) 

x H 2 0 -> x H 2 + 0.5x 0 2 

(30) 
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Figure 12. — Diagram of the carbothermal reduction process with identified features. 


Orbital Technologies Corporation (ORBITEC) has developed a module to extract oxygen using the 
Carbothermal Reduction method discussed above. The modeled system is shown in Figure 12. The steps in the 
process are as follows: 

• Scoop the regolith from the surface and transport it to the module. 

• Separate the large particles from the regolith using a vibrating screen. 

• Transport the regolith to the regolith hopper then to the reactor chamber using two augers. 

• Heat the regolith to 1900 °C by using concentrated solar energy, and circulate methane. 

• Extract the processed regolith and transport the produced carbon monoxide and hydrogen gas to the 
methanation reactor. 

• Mix additional hydrogen with the products in the methanation reactor to produce water and methane. 

• Methane is returned to the system. 

• Electrolyze the water to produce hydrogen and oxygen. 

• The hydrogen is returned to the system and the oxygen is stored. 

The first step in the process is to gather and transport the regolith to the module. This analysis considers a belt 
scoop to accomplish both tasks. The regolith would then be screened to remove large particles. This would ensure 
that the system did not jam. The next step involved transporting the regolith to the reactor chamber, which was done 
by two augers. Once in the reactor chamber, concentrated solar energy would be used to melt the regolith. The 
chamber could process 7 melts simultaneously. After the melt(s) reach the required temperature, methane gas would 
be introduced and the first reaction would take place (Eq. (28)). The processed regolith could then be scooped out 
using a rake, while the reactant gasses (CO and H 2 ) would enter into the methanation reactor (Eq. (29)). Water 
would then be extracted and cleaned of any impurities. Finally, oxygen is extracted from the water by an 
electrolyzer. 

To determine the total power requirement of the entire system, the power requirement of each individual 
component was either calculated or estimated. Since it was known that the desired oxygen production is 
1000 kg/year and that the typical oxygen yield is 15 to 20 wt%, the minimum mass flow rate of the regolith can be 
determined (Ref. 10). It was assumed that the lunar location could be at any site, so therefore the carbothermal 
module would only be operational approximately half of the year (4380 hr). The mass flow rate of the regolith 

( M r ) can be determined from Equation (31): 


M 


r 


mg 

t s (yo ) 


(31) 


Where m 0 is the desired oxygen production per year (kg), t s is the total amount of time of available sunlight 
(hours), and y 0 is the oxygen yield from the regolith. The mass flow rate of the regolith was used to determine the 
amount of power required at each stage. The first stage where power is needed is to scoop and lift the regolith. To do 
this, a belt scoop is utilized. The power of the scoop came from lifting the soil as well as the drag resistance against 
the regolith. It is assumed that the scoop is in contact with the reqolith for 1 m. The distance to lift the regolith was 
determined by the height of the carbothermal module plus its distance above the surface. The electric motor 
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efficiency was assumed to be 90 percent. The power requirement for scooping the regolith (P ra ) and lifting the 
regolith (P w ) can be determined by Equations (7) and (8). 

Several of the terms are identical to the hydrogen reduction case: r\ em is the electrical efficiency of the motor, \i r 
was the coefficient of friction of the regolith (assumed to be similar to sand or 0.4), and d rs was the distance the 
scoop came in contact with the regolith (assumed 1 m). The carbothermal reactor the height (h r ) was 1.37 m) 

(Ref. 10) and h rs was the height of the extending belt scoop (assumed 1 m). The second step was to separate out the 
large particles of regolith by using a vibrating screen. The largest particle size (</*) was assumed to be identical to the 
particle size required for the hydrogen reduction process, which was determined to be 0.84 mm. The power required 
to operate the vibrating screen can be determined by Equation (12). 

After being screened, the regolith would then be processed. Inside the system, several augers and actuators were 
used to move the machine components and/or the regolith. 50 W were assumed as the power requirements for the 
machine movements. This should provide some margin based on the current power requirement estimates. The 
ORBITEC module produced regolith melts that were approximately 7.5 cm across and were roughly hemi-spherical 
(Ref. 10). The number of melts per hour could be determined based on the volume of each melt, the flow rate of the 
regolith, and the density of the regolith, which was approximately 1000 kg/m 3 . This number was calculated as 
13.8 per hour, and dictates how much time could be allotted for each melt to maintain the yearly oxygen production 
quota. The thermal power needed to heat the regolith to the desired temperature ( Q h ) could be found by 
Equation (32): 


Where n m was the number of melts being processed in the module, p r was the density of the regolith, c p was the 
specific heat (1646.8 J/kg*K based on the desired temperature for the reaction), v w was the volume of a single melt, 
T m was the temperature of the melt, T r was the temperature of the regolith, and t m was the time allowed for the 
regolith to reach the desired temperature. The desired temperature of the melt was 1900 °C (Ref. 10). The 
temperature of the regolith was the same as the temperature on the surface of the moon during the day or 384 K 
(Refs. 1 and 7). The relation for the specific heat of the regolith came from a curve fit of experimental data on the 
regolith (Ref. 5). This relationship is described by Equation (1). 

The regolith melts are heated in a bed of regolith that utilizes the surrounding regolith as insulation (Fig. 13). 
Due to this thick layer of insulation, it was assumed that no heat was lost through the regolith. This structure limits 
radiation losses to the surface of the regolith. 

The equation for radiation heat transfer was utilized along with the number of regolith melts being processed. 
Since the processed regolith needed to be solidified before it could be extracted, it was assumed that the melt would 
be held at the processing temperature half of the time and allowed to solidify the other half. Equation (33) was used 
to calculate heat losses through radiation (g,): 



(32) 



(33) 


r~ Melted 

r~ Cool Regolith / Re 9 olith j— insulating Regolith 



Figure 13. — Diagram of Regolith melt insulated by surrounding 
Regolith. 
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Where A s was the surface area exposed to radiation loses, ct was the Stefan-Boltzman constant with a value of 
5.67x 10 8 W/m 2 K, T s was the temperature of the surrounding (assumed to be 270 K), e c was the emissivity of the 
container of the reactor (assumed to be chromium with a value of 0.14), e* was the emissivity of the insulation 
material around the reactor (assumed to be nickel with a value of 0.17), and rii was the number of insulation layers 
(assumed 50). A well-insulated system would only require an estimated 1 or 2 W to keep the system at daytime 
temperatures. Other than radiation loses, the regolith would also absorb heat due to the heat of reaction of the 
reactants ( Q r ). This was determined by summing the heats of reaction for each of the substances depending on the 
percentage of each of the substances in the regolith (discussed further in the next section). The next step in the 
process was to have the reactant gasses enter the methanation reactor (combining CO and H 2 to form CH 4 and H 2 0) 
as well as clean any impurities contained within the gases. It was assumed the energy produced in the methanation 
reaction met the energy requirement for the clean-up modules. The final power requirement was the electrolyzer ( P e ) 
for the H 2 0 to produce H 2 and 0 2 . This was calculated by using the 3520 W/kg/h theoretical power needed to 
separate the components and the efficiency of the electrolyzer, which was assumed to be 0.72, (Ref. 7) see 
Equation (23). The total thermal and electrical power requirements were determined by summing the individual 
component power requirements discussed above. 


Results of Sizing for Carbothermal Reduction 

Based on the desired 1000 kg/year oxygen production rate with an assumed 15 percent yield, the mass flow rate of 
the regolith was determined to be approximately 1 .52 kg/h. This flow rate was used to determine the power 
requirements of each component of the system based on one melt being processed at a time. The thermal power 
requirements can be found in Table 4, and the electrical power requirements can be found in Table 5. A visual 
representation of the breakdown of power requirements (both thermal and electrical) can be seen in Figures 14 and 15. 


TABLE 4.— BREAKDOWN OF THERMAL POWER REQUIREMENTS 


System 

Power, 

W 

Radiation heat loss 

272.5 

Power to heat the regolith 

2491.2 

Total heat of reaction 

5555.5 

Total thermal power 

8319.2 


TABLE 5.— BREAKDOWN OF ELECTRICAL POWER REQUIREMENTS 


System 

Power, 

W 

Power to lift the regolith 

0.0011 

Power to scoop the regolith 

0.0003 

Power of the vibrating screen 

0.3651 

Power of the actuators and augers 

50.00 

Power of the electrolyzer 

1754.1 

Total electrical power 

1804.4 
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■ Heat loss to surroundings 

■ Regolith heating 
Heat of reaction 



Figure 14. — Breakdown of thermal power requirements for 
carbothermal reaction. 


Operate electrolyzer 

■ Process Regolith 

(lift, scoop, vibrate, separate) 



Figure 15. — Breakdown of electrical power requirements for 
carbothermal reaction. 
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To determine the total heat of reaction, each compound in Table 1 was analyzed for their individual heats of 
reaction with CH4. These reactions were assumed to be at room temperature, since the difference equated to less than 
0.3 percent for Si0 2 reacting at 1900 and 25 °C (room temperature). The heats of reaction were calculated to be 
approximately 841 and 839 kJ/mol Si0 2 , respectively. The calculated heats of reaction for each of the major substances 
are shown in Table 6. 

The total heat of reaction depended on the amount of each substance (in Table 6) that was present in the regolith 
(the percent composition can be found in Table 1). The total estimated heats of reaction for the Highland and Mare 
regions were 5,556 and 5,135 W, respectively. To ensure that the calculations did not bias one area over the other, 
5,556 W was used for the total heat of reaction. 

The calculated power requirements above were based on the ORBITEC carbothermal regolith reduction module 
running one batch (melt) at a time. However, the module could run 7 melts (batches) simultaneously (Ref. 10). Since 
each melt had a specific size and the mass flow rate of the regolith needed to meet a specific value, the time allowed 
for each melt process could not exceed a specific value. It should also be noted that under certain cases the required 
process time could be unrealistic. The maximum time allowed for each melt under different numbers of 
simultaneous processes can be seen in Figure 16. 

It was assumed that half the total amount of time would be spent on heating the regolith. However if more time 
was required, more modules could be simultaneously operated in order to maintain the desired production rate of 0 2 . 
Also due to the additional melts, more thermal power would be required. Based on the increased thermal power for 
the additional melts being processed, the total power required for various numbers of melts can be seen in Figure 17. 


TABLE 6 .— HEAT OF REACTION FOR METAL OXIDE 
WITH CH4 TO PRODUCE CO AND H 2 


Reactant 

Heats of reaction, 
kJ/mol 

Si0 2 

838.84 

AI 2 O 3 

1567.91 

FeO 

236.07 

CaO 

598.99 

MgO 

565.67 

Ti0 2 

872.14 




Figure 16. — Maximum time available to process each 

batch-dependent upon the amount of oxygen produced Figure 17. — Total thermal power requirements based on 

per year (1 000 kg/year) and utilizing only one module. number of simultaneous batches using one module. 
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The power consumption for the Carbothermal Reduction can be compared to the power consumption of the 
Hydrogen Reduction method. Using the scenario of only 1 melt operating at a time, the total power requirement for 
the carbothermal process came to 10,124 W. With a solar intensity of 1307 W/m 2 (approximated solar intensity on 
July 9th at 180°), a solar concentrator with an area of 7.81 m 2 would be needed for a single melt; an area of 9.35 m 2 
would be needed for 7 melts, with radii of 3.15 and 3.45 m, respectively. For the carbothermal system it was 
assumed that the electrical power needed to run the system would be provided by solar arrays. Table 7 shows a 
representative sampling of terrestrial rigid panel solar arrays along with the total area and mass needed for the arrays 
(based on the power output of each array and the total electrical power needed) (Ref. 11). The power consumption of 
the Hydrogen Reduction method was approximately 7,400 W. However, much of that power was used for the 
Stirling engine. Removing the heat required for the engine left a power requirement of approximately 3,71 8 W based 
on an equivalent 0 2 mass flow rate. 


TABLE 7,— REPRESENTATIVE TERRESTRIAL RIGID PANEL SOLAR ARRAY AREA AND MASS (REF. 11) 



Mass, 

Area, 

Power, 

Number 

Total area, 

Total mass, 


kg 

m 2 

W 

needed 

m 2 

kg 

M55 

5.5 

0.43 

55 

32.8 

14.0 

181.9 

SX-170B 

15.0 

1.26 

170 

10.6 

13.4 

159.2 

ES-190-SL 

18.2 

1.50 

190 

9.5 

14.2 

173.1 

GEPV-200 

17.7 

1.44 

200 

9.0 

13.0 

160.0 

1KA-GSA 

13.7 

0.95 

60 

30.1 

28.6 

412.0 

KC130GT 

15.0 

0.93 

130 

13.9 

12.9 

208.2 

PV-MF170EB4 

15.5 

1.26 

170 

10.6 

13.4 

164.5 

PW1650-165 

18.0 

1.36 

165 

10.9 

14.9 

196.8 

HIP-205BA3 

14.0 

1.18 

205 

8.8 

10.4 

123.2 

S-170-SPU 

15.5 

1.26 

170 

10.6 

13.4 

164.5 

ND-208U1 

21.0 

1.63 

208 

8.7 

14.1 

182.2 

175-PC 

17.3 

1.32 

175 

10.3 

13.6 

178.1 

SW165 

15.0 

1.30 

165 

10.9 

14.3 

164.0 

STP170S 

15.5 

1.28 

170 

10.6 

13.6 

164.5 

SW160 

20.0 

1.30 

160 

11.3 

14.7 

225.6 

Average 

15.8 

1.23 

160 

13.2 

14.6 

190.5 


Conclusions 

The results of the baseline system produced just under 0.6 kg of oxygen per hour with a mirror size of less than 
5 m diameter. For a full 14-day mission, where the production plant is in continuous sunlight, approximately 195 kg 
of oxygen would be produced. This translates into a yearly production rate of approximately 2300 kg of oxygen. The 
system size should also enable it to be packaged into a single launch vehicle providing a system configuration that 
can be autonomously deployed and operated. If greater oxygen production rates are required the system can be 
scaled up, requiring a larger mirror size for operation. Another alternative is to consider these oxygen production 
systems as modular. If an increased rate of production is needed then two or more systems can be utilized. 

As more detail is incorporated into the various components that make up the system, the sizing and oxygen 
production rates will likely change. However, the results produced by this analysis should be representative of the 
systems’ potential. The analysis done demonstrated that more power would be required to produce 1000 kg/year of 
oxygen using the carbothermal reduction method as compared to the hydrogen reduction method. However, there is 
the advantage of utilizing all of the regolith instead of just the ilmenite that is used for hydrogen reduction. 
Furthermore, a solar concentrator with a diameter less than 4 m could provide the thermal power needed. The 
electrical power (1804 W) would be provided from solar panels. 

An integrated solar dynamic/Stirling heat engine system or a hybrid solar/photovoltaic are excellent candidates 
for powering a lunar oxygen production facility utilizing respectively, either hydrogen reduction or carbothermal 
reduction technologies. The flight system should provide for an efficient, compact design that can be designed as an 
integrated stand-alone unit. This would enable the system to be utilized for a number of different applications from a 
main base to a remote camp. It also can provide a basis for a standardized, modular oxygen production system 
design. 
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photovoltaics are compared for the production of electricity. The electricity produced is utilized to operate the equipment needed in the 
oxygen production process. The initial oxygen production method utilized in the analysis is hydrogen reduction of ilmenite. Utilizing this 
method of oxygen production a baseline system design was produced. This baseline system had an oxygen production rate of 0.6 kg/hr with 
a concentrator mirror size of 5 m. Variations were performed on the baseline design to show how changes in the system size and process 
(rate) affected the oxygen production rate. An evaluation of the power requirements for a carbothermal lunar regolith reduction reactor has 
also been conducted. The reactor had a total power requirement between 8,320 to 9,961 W when producing 1000 kg/year of oxygen. The 
solar concentrator used to provide the thermal power (over 82 percent of the total energy requirement) would have a diameter of less than 4 
m. 
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